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Rejoinder on: Thermostatistics of Overdamped 
Motion of Interacting Particles 

In their Reply [H to our Comment Andrade et al. 
state that we have "chosen to categorically dismiss their 
elaborate and solid conceptual approach without employ- 
ing any concepts or tools from Statistical Mechanics". In 
this Rejoinder we will show that contrary to the state- 
ment above, approach employed by Andrade et al. is 
neither "solid" nor "sound". Unfortunately, because of 
the one page restriction imposed by PRL we could not 
address all of the flaws of the original paper. Therefore, 
we are grateful to Andrade et al., for giving us an oppor- 
tunity to further elaborate on our Comment. Indeed, in 
our Comment we have shown that at T = the model 
studied by Andrade et al. Q can be solved exactly. In 
their Reply Andrade et al. argued that the thermody- 
namic limit employed by us to obtain the exact solution 
was "less physical" than their simulation in which 800 
vertices were confined in a strip with periodic boundary 
condition in one direction and an external parabolic po- 
tential in the other. To what extent this is more realistic 
than our thermodynamic limit calculation is very ques- 
tionable, but lets leave it at that. 

1. Let us first briefly address the point of thermody- 
namic limit. Andrade et al. surely realize that sta- 
tistical mechanics, strictly speaking, is valid only in 
the thermodynamic limit. Without this limit there 
is no equivalence between different ensembles. For 
particles confined inside a parabolic potential well 
there are two ways that one can perform the ther- 
modynamic limit (a) one can scale the vertex charge 
q so that it is proportional 1 / \/N or alternatively 
(b) one can leave q=l and scale the confining po- 
tential so that its stiffness a (x N. In either case the 
reduced one-particle distribution function will then 
converge to the same function in the limit, N ^ oo. 
The way that Andrade et al. performed their sim- 
ulations, there is no thermodynamic limit and for 
each different N the distribution will be different. 
In thermodynamic limit, the model of Andrade et 
al. can be solved exactly for any T. In the Com- 
ment, we only presented the solution for T = 0, but 
it can easily be extended to any T, showing that the 
system of Andrade et al. is always described by the 
usual Boltzmann-Gibbs statistical mechanics. 

2. But lets us for a moment forget about the ther- 
modynamic limit. In this case, mean-field is no 
longer exact and the particle correlations must be 
taken into account. We first note that a finite num- 
ber of particles confined in a periodic potential at 
T = will crystallize - of course, because of the 
confinement, the crystal will have many defects. 
It is, therefore, meaningless to assign a uniform 
functional form to density, which in reality will be 
made of a sum of delta-functions. It is like saying 
that NaCl crystal is described by a constant den- 
sity, with no Bragg peaks or anything related to the 



crystal structure. At T = the only way to make 
sense of a differential equation like eq.(l) of Reply 
is if iV — > oo. 

But even allowing the incorrect coarse-graining pro- 
cedure of Andrade et al., we still see that the 
density distribution will not be given by a q- 
exponential. A q-exponential is a smooth function, 
so that the density must go continuously to zero. 
This is not the case for the system of Andrade et al. 
In our Comment we saw that the particle density 
drops discontinuously to zero. For finite N, the par- 
ticle correlations will change slightly the form of the 
distribution function from the one obtained in the 
thermodynamic limit, but they will not change the 
fact that there is a discontinuity in density. This 
can be clearly seen from the simulation of Andrade 
et al. presented in Fig.l of their Reply. The figure 
clearly shows a discontinuous jump in density. Be- 
cause of the way that Andrade et al. performed the 
binning to construct their histogram, the jump has 
been smoothed out somewhat. We have rerun the 
simulations of Andrade et al. for exactly the same 
parameters and have plotted the data in Fig. 1. 
The figure clearly shows a discontinuity in density, 
after the leftmost and the rightmost points of the 
histogram the density is identically zero. This dis- 
continuity is predicted by our theory, and is absent 
within the Tsallis formalism. We conclude again 
that the particle distribution at T = has nothing 
to do with Tsallis statistics. 
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FIG. 1: The density distribution at T = for iV = 800 inside 
tlie parabolic potential well with a = 10~^/o. The simulation 
is identical to the one done by Andrade et al. After the two 
extreme points on the left and on the right, the density is 
identically zero. 



4. But again let us give a benefit of the doubt to 
Andrade et al. and ask the authors the follow- 
ing question. For T = it is very simple to find 
the position of mechanical equilibrium for N—1 or 
N=2,3,4,5...10 etc. We can also calculate the par- 
ticle distribution in the thermodynamic limit, as 
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was done in our Comment. In all these cases, the 
equilibrium is simply determined by the force bal- 
ance. No statistics is involved. Our question is this: 
for exactly what value of N do authors expect that 
Tsallis entropy will start determining the particle 
distribution at T=0? 

5. Now suppose that we raise the temperature. At 
some point the crystal will melt. However, the fluid 
state will still have strong correlations between in- 
dividual particles. There is more than a hundred 
years of history, going all the way back to the pi- 
oneering works of Debye, Ornstein, Zernike, and 
Kirkwood on how the particle correlations can be 
calculated. There are thousands of liquid-state the- 
orists who have dedicated their lives to understand- 
ing statistical mechanics of liquids — all this work 
simply ignored by Andrade et ai. In fact the corre- 
lations and their effects on the density distribution 
can be calculated in the framework of the usual 
Boltzmann-Gibbs statistical mechanics. There are 
theories, such as the Hypernetted Chain Equation, 
Roger- Young Equation, SCOZA, the Density Func- 
tional Theory, etc. which were derived to do exactly 
this. There is no reason to introduce any fitting 
parameters through Tsallis entropy. This does not 
lead to any new physical understanding, only to 
curve fitting. In fact, in their paper Andrade et al. 
did not provide any reason or indication why they 
believe that the standard Boltzmann-Gibbs statisti- 
cal mechanics will not apply to the system of ver- 
tices studied by them. 

6. In their Reply, Andrade et al. say that we did 
not question their non-linear Fokker-Planck equa- 
tion (eq. 1 of Reply) on which are based all of 
their arguments. Indeed, because of the one page 
restriction we could not get into detailed discus- 
sion of the original paper, concentrating only on 
our exact solution for T = 0. In fact, we do ques- 
tion the use of this non-linear Fokker-Planck equa- 
tion to study the dynamics of interacting particles. 
This is an incorrect equation to use for this system. 
It is well known that the dynamics of interacting 
particles is governed by the self-consistent Nernst- 
Planck equation. Again, Andrade et al. got lucky 
because for the particles interacting through a mod- 
ified Bessel function, the stationary solution (after 
adjusted with the fitting parameter as was done by 
Andrade et al.) is quite close to the correct solution 
of Nernst-Planck Equation. Thus, the asymptotic 
dynamics derived using the incorrect Fokker-Planck 
equation actually looks quite reasonable. In fact it 
is simple to see that eq. (1) used by Andrade et al. 
is incorrect. Using it Andrade et al. derived the 
unlucky eq. 13 of their PRL, from which all the 
incorrect discussion followed. From this equation 
Andrade et al. concluded that the stationary state 
for particles interacting by any short-range force is 



always parabolic - a Tsallis q-Gaussian of entropic 
index v = 2. From our exact calculation, it is clear 
that the parabolic form is very special for the inter- 
action potential of the modified Bessel form. Other 
potentials will not lead to the parabolic form. We 
should recall that Andrade et al. claim that their 
theory applies to particles interacting by an arbi- 
trary short-range potential. As a simple demon- 
stration that this is not true, we have simulated a 
Id systems of particles interacting through a short 
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range potential V{x) = \e~^ . In Fig. 2, we see 
that once the system has relaxed to equilibrium, 
the distribution is not the inverse parabola, but a 
much more complicated function. One should also 
note the characteristic discontinuity in density, con- 
tradictory to Tsallis statistics. We conclude again 
that the non-linear Fokker-Planck equation of An- 
drade et al. can neither account for the stationary 
state nor for the dynamics of approach to equilib- 
rium. 




0.1- 

0.0 ^ \ \ \ \ \ 

-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

X 



FIG. 2: The density distribution at T = for particles in- 
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teracting by a short range potential V(x) ~ ^ confined 
in a parabolic trap. To get better statistics and faster run 
time, we studied a Id model, with A'^ — 1000, a = 0.05. Note, 
however, that because of the periodic boundary condition the 
model of Andrade et al. is also effectively Id. Also note the 
discontinuous drop in density. 



7. Andrade et al. have misinterpreted the interparti- 
cle correlations present in a finite system for exis- 
tence of a new kind of entropy. The fact that the 
system studied by Andrade et al. is described by 
the usual Boltzmann-Gibbs statistical mechanics, 
follows simply from the study of velocity distribu- 
tion in equilibrium. In Fig. 3 we show a result of a 
microcanonical (constant energy) simulation. The 
system starts from an arbitrary initial condition 
and evolves in accordance with Newton's equations 
of motion. After some time it relaxes to thermody- 
namic equilibrium. In Fig. 2 we plot the histogram 
of particle velocities together with the 2d Maxwell- 
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Boltzmann distribution, 
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We see a perfect agreement between the 
Boltzmann-Gibbs statistical mechanics and 
the simulations, without any adjustable parame- 
ters. In this example, the equilibrium temperature, 
T = {v'^) /2, is very low T = 0.3 and Andrade et al. 
claim that the system should be described by the 
Tsallis statistics. In fact we do not see any trace of 
q-exponentials at this or at any other temperature. 

In their Reply Andrade et al. state that "neither 
the Maxwell-Boltzmann nor the Tsallis thermostat- 
ics are contrary to classical mechanics theory". We 
respectfully disagree with this statement. Indeed 
the Boltzmann-Gibbs statistical mechanics is in full 
agreement with the classical mechanics. However, 
there is not a single classical system of particles 
interacting by a short-range potential (for long- 
range forces see [1] and the references therein) that 
evolves to "Tsallis" equilibrium. The model stud- 
ied by Andrade et al. is precisely the case in point. 
Therefore, we stand by our original statement: "the 
density distribution of particles in contact with a 
reservoir at T = has nothing to do with the Tsallis 
statistics, and everything to do with the Newton's 
Second Law." 



FIG. 3: Velocity distribution. Points are from molecular- 
dynamics simulation, for exactly the system studied by An- 
drade et al, N = 800, Ly = 20, a = 10"^/o. Solid curve 
is the Maxwell-Boltzmann velocity distribution. There are 
no adjustable parameters. The temperature is T = 0.3. At 
this temperature, Andrade et al. claim that the system is de- 
scribed by the Tsallis statistics. As is clear from the graph, 
we see no trace of Tsallis q-exponentials in the velocity dis- 
tribution function. 
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